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ABSTRACT 

The purpose of this paper is to present the resuhs of a theoretical study 
pertaining to the feasibihty of PIAA units using Deformable Mirrors. We begin by 
reviewing the general derivation of the design equations driving PIAA. We then 
show how to solve these equations for square apertures and show the performance 
of pure PIAA systems in the ray optics regime. We tie these design equations 
into the study of edge diffraction effects, and provide a general expression for the 
field after a full propagation through a PIAA coronagraph. Third, we illustrate 
how a combination of pre and post apodisers yields to a contrast of 10^^° even 
in the presence of diffr active effects, for configuration with neither wavefront 
errors or wavefront control. Finally we present novel PIAA configurations over 
square apertures which circumvent the constraints on the manufacturing of PIAA 
optics by inducing the apodisation with two square Deformable Mirrors (DM). 
Such solutions rely on pupil size smaller than currently envisioned static PIAA 
solutions and thus require aggressive pre and post-apodizing screens in order 
to mitigate for diffractive effect between the two mirrors. As a result they are 
associated to significant loss in performance, throughput in particular. 

Subject headings: instrumentation: high angular resolution 
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Introduction 



Direct detection and spectral characterization of exo-planets is one of the most exciting 
challenges of modern astronomy. Over the past decade, the field of high contrast imaging 
has been extremely active, and a variety of solutions for space based imaging of earth twins 
have been proposed and tested. Among thes e concepts, th e Phase Induced Amplitude 



Apodisation coronagraph, first introduced by 



GuyonI (120031 ) . is a promising solution since 



it makes most of the photons collected by the primary mirror of the telescope available for 
planet detection and characterization. This technique is based on two aspherical mirrors 
that redistribute the light in the pupil plane of a telescope so that it follows a given 
amplitude profile, creating a Point Spread Function (PSF) with contrast levels close to 
10^^*^. Because all the light collected is remapped using these mirrors, such a coronagraph 
has virtually no throughput loss. As a consequence, the angular resolution is undiminished 
and is close to the diffraction limit of IX/ D, a feature comparable to the performance s 



of a p hase mask coronagraph, such as the Optical Vortex Coronagraph fiMawet et al. 



( 120101 )). but without any transm issive optics. Unfortunate 



on a geometries 



optics analysis (ITraub fc Vanderbei 



(120031) 



y, these conclusion s are 



3ased 



Martinache et al. 



(120061 ) ). 



Vanderbei! ( 120061 ) discovered that when a full diffraction analysis is applied to the system, 
the envisioned designs do not yield the expected level of extinction, but rather, at best, 
10~^ contrast. 



In order to mitigate these diffractive effects, iPluzhnik et al.l (120051 ) devised a method 
using pre -and post-apodisers to retrieve a 10~^° contrast even while accounting for 
diffraction. They verified, using a very accurate numerical integrator, that apodised PIAA 
coronagraphs could be designed to have performances close to the ones predicted by ray 
optics. In the absence of apodizing screens, another avenue to mitigate for propagation 



induced artifacts is to treat them as wavefront errors and compensate for them with one 
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or two Deformable Mirrors (iTrauger fc Traubl (120071 ). iPueyo et al.l (l2009a( ) jGuyon et al. 
(120101)) • The first solution, apodizing screens, will reduce the throughput of the coronagraph, 
and thus the efficiency of the instrument, but preserve the spectral bandwidth. The second 
alternative, wavefront control, will preserve throughput but lead to a reduced spectral 
bandwidth since the diffractive oscillations are highly chromatic and cannot be fully 
compensated by deformable mirrors. In this communication we treat only the former design 
approach, where diffractive effects are mitigated by pre- and post-apodizers. 



In this pa p er we extend the theoretical framework of 



Pluzhnik et al 



Vanderbeil (120061 ) and 



(|2005[ ) by developing an analytical formalism for the diffractive effects on 



arbitrary geometries. These results are not only useful for a qualitative understanding 



of pupil mapping coronagrap^ 



Balasubramanian et al 



aberrations. 



Pueyo et al. 



control systems. 



IS, but are also critical for the design of such solutions. 



( 2010 ). for numerical studies pertaining to their sensitivities to 



1120 



Krist et al. 



^and for understanding how to integrate with wavefront 



( 120101 ). We apply these findings to pursue alternate designs 
of PIAA coronagraphs, and focus on the particularly interesting case of square apertures. 
Indeed, if the apodisa tion function is written as th e product of two separable, one- 



dimensional functions, 



Nisenson fc Papalioliosl ( 120011 ). very deep contrast can be achieved 



on the four diagonal directions of the Point Spread Function (PSF) with only mildly 
deformed PIAA mirrors. This configuration implies that if the telescope aperture is circular, 
the square sub-aperture needs to be inscribed in the telescope pupil, leading therefore to a 
throughput loss of a factor of 2/7r. Since any exo-planet dedicated coronagraph will have to 
be equipped with two sequential DMs to reduce the stella r halo due to i mperfe ction in the 



optics at mid-spatial frequencies (IShaklan fc Greenl ( 120061 ): 



Pueyo et al. 



(l2009al )). and since 



current state of the art pi ezo-electric or MEMS (M i croelectromechanica. 



all have a square formats ( iTrauger &: Traubl ( 120071 ) 



Give 'on et al. 



Systems) devices 



mom , it is sensible to 



study square apertures inscribed in a circular telescope aperture. Alternatively, one could 
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envision a design for which the circular pupil is inscribed in a square deformable mirror. 
This would increase the net throughput of the system, and the Inner Working Angle (IWA) 
would become constant for all azimuthal angles of the PSF. However, the Outer Working 
Angle (OWA) would not take advantage of all the actuators on the Deformable Mirror in 
the diagonal direction and the PIAA deformations would become larger. Adressing such a 
paradigm is beyond the scope of this paper, and we focus in this communication on the 
case where the square pupil is inscribed in the circular aperture, t he cas e that is wide ly 



i mpleraented in current hig h contrast testbeds (ITrauger &: Traubl (120071 ) 



(|2009a[ ) 



Guyon et al 



Pueyo et al 



torn ). 



We begin by reviewing the general derivation of the defining PIAA equations and 
show how to solve these equations in the simplified case of square separable apertures. 
We present the instrument response in that case, and study the performance of this pure 
PIAA systems under the ray optics approximation. In a second section we revisit the 
analysis considering edge diffraction effects and provide a general expression for the field 
after a full propagation through any arbitrary PIAA coronagraph. In the remainder of the 
paper we apply these general results to the design and trade-offs associated with PIAA 
coronagraphs over square separable apertures. We present numerical results that use our 
new diffractive approximation and illustrate how a combination of pre- and post-apodisers 
restores contrast below 10^^° even in the presence of diffractive effects. Finally, using the 
result of this theoretical study, we present novel PIAA designs over square apertures which 
can potentially relax the constraints on the manufacturing of PIAA optics or completely 
circumvent manufacturing by inducing the apodisation with two square Deformable Mirrors 
(DM). 
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2. Design of Pupil remapping systems 



2.1. General equations 



The Phase Induced Amplitude Apodization (PIAA) coronagraph was fist introduced 
by iGuyonI (120031 ) as an optical apparatus composed of two refiective surfaces that remapped 
the distribution of the light in the pupil of a telescope into an apodisation function which 



yields a high contrast PSF. 



GuyonI (120031 ) showed that the mirror surfaces co uld be derived 



Traub &: Vanderbeil (120031 ) subsequently 



as the solution of a single differential equation and 
developed this result showing that the relationship between prescribed apodisation and 
mirrors surfaces could also be written as a set of two coupled partial differential equations. 
Here we rederive these ray-optics based design equations for arbitrary geometries and 
illustrate how to solve them using square apertures. 

Consider the arrangement and coordinate systems shown in Fig. [H Note that in this 
paper we only consider the on-axis pupil-to-pupil configurations leaving aside the off-axis 
focal-to-focal case, whose surfaces cannot be derived by direct integration of the systems of 
partial differential equations. For now we do not assume anything about the geometry of 
the two remapping mirrors: Ml and M2 are two domains of of arbitrary shapes. We 
introduce the independent variables (x, y) that span the domain of Ml and (x, y) that span 
the domain of M2. The height of Ml is thus given by the function h{x, y) and the height of 
M2 is given by h{x,y). 

Given the distance between the two mirrors, Z, the optical path length, Q, between 
any two points (x, y) and (x, y) is given by 



Z Z 

Q{x,y,x,y) = h{x,y) - — + S{x,y,x,y) + — - h{x,y) 

= S{x,y,x,y) + h{x,y) -h{x,y) 
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where S is the free space propagation within the apparatus 



S{x, y, X, y) = \/{x- + {y - y f + {h{x, y) - h{x, y)y. (2) 



Following a two-dimensional generalization of the approach in IVanderbei fc Traub 



( l2005l ). the pupil mapping solution for the mirror heights is a one-to-one mapping between 
points in the domain of Ml to points in the domain of M2. For instance, if we consider the 
coordinates {x, y) on Ml as independent variables, then the pupil mapping solution maps 
rays to specific points {xo,yo) via two functions 

Xo = fiix,y) (3) 
yo = f2ix,y). (4) 

Alternatively, since the pupil mapping solution must be reversible, we can consider the 
points (x, y) as independent variables and find the inverse mapping, 

Xi = gi{,x,y) (5) 
yi = g2{x,y) (6) 

where gi and g2 are the inverse functions of /i and /2. Typically we design the pupil 
mapping system by selecting the inverse functions gi{x,y) and g2{x,y) since our goal is a 
specific amplitude profile at M2. 

We design a pupil remapping unit via ray optics by choosing the heights of each mirror, 
h{x,y) and h{x,y), so that the output amplitude profile follows the desired apodisation 
A{x,y). This requires that the following three constraints be satisfied: 



1. The intensity profile after M2 matches the square of the desired apodisation A{x, y)^. 
The relationship between this profile and the mirror surfaces has to satisfy local 
energy conservation between Ml and M2. 
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2. For each ray between {xi,yi) and {xo,yo), Fermat's principle must be satisfied. 

3. When the apparatus is illuminated by uniform on-axis light, the outgoing wavefront is 
flat. This means that, in the absence of diffractive effects, the phase at M2 is constant 
when the wavefront hitting Ml is flat. 

In the remainder of this section we show how, for any arbitrary geometry, these constraints 
translate into formal conditions on the heights of each mirror under a ray optics approxi- 
mation. The set of equations obtained defines the design of a PIAA unit. 



2.1.1. Energy remapping 

Constraint [1] is simply a statement of energy conservation. This translates directly into 
an equation for the remapping functions gi and g2 (or, equivalently, /i and /2). Fig. |2] 
illustrates that the energy in a small square surface element of Ml must equal the energy in 
its remapped image at M2. Because of the aspherical nature of the PIAA mirrors the area 
of the each element varies as a function of location at M2 and thus, in order to conserve 
energy, the light intensity in each square changes accordingly. More formally, this implies 
that the energy contained in a differential area dxdy at Ml must remain constant when 
traveling through the set of aspherical PIAA optics. Considering (£, y) as independent 
variables, we can view the inverse mapping as a coordinate transformation, given by Eqs. 
and [6l allowing us to use the differential area formula, 

dxdy = det( J (x, y)) dxdy (7) 

where J{x,y) is the Jacobian matrix of the change of coordinates {x,y) to {xi,yi). Since 
the remapped rays correspond to an amplitude profile A{x,y), the energy in a bundle of 
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rays at M2 is given by A{x,yy. Consequently, the energy remapping constraint in Eq. [7] 
can be written as 

dgidg2 dgidg2 

A{x,y) 



dx dy dy dx 

where, again, gi{x,y) and g2{x,y) are the remapping functions in Eqs. [5] and |6l 

The desired apodisation thus does not directly lead to the mirror surface heights but 
only to the location of the incident rays. This implies, as we will study in the next section, 
that the output field strongly depends on how well the ray optics approximation holds. 
Eq. [8] is a two-dimensional version of th e well known n on-linear partial differential equation 



known as the Monge-Ampere equation. 



Mongd (jl78ll ). While there is no known general 



solution to this equation for gi and g2, we show that, by using a square aperture we can 
find solutions if we restrict the apodisation to be the tensor product of two one- dimensional 
functions. 



2.1.2. Mirror Heights 



To find the heights of Ml and M2 we follow the approach in lVanderbei &: Traubl (120051 ) 
and satisfy Fermat's principle (constraint [2]) by minimizing the optical path length of each 
ray from to (xo,yo)- After some algebraic manipulations this provides the equation 

governing the height of Ml at each (x, y) across the mirror. 



fi{x,y) - X 



dy 



x,y)) 



S{x,y J1J2) + {h{x,y) - h{f 1J2)) 

f2{x,y)-y 

S{x,y J1J2) + {Kx,y) - h{fij2)) 



(9) 
(10) 



For the height of M2 we use constraint [31 whic h states that the phase of ;he wavefront 
be constant across the exit pupil. It was shown in iVanderbei fc Traubl (120051 ) that this 
translates into the condition that each ray between the pupil remapped points {xi,yi) and 
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[xoiUo) must have constant optical path length. From Eq. [T]this means that 



= S{x,y J1J2) + (h{x,y) -h{f I J2) 
= 2Z. 



'IV. 



Taking the partial derivative of the optical path length with respect to (x, y) and 
setting it equal to zero leads to expressions for the height of M2, 

X- giix,y) 



^ix,y) 



Sigi,g2,x,y) + {h{gi,g2) - h{x,y)) 

y -92{x,y) 

Sigi, g2, X, y) + {h{gi, g2) - h{x, y)) 



(12) 
(13) 



The resulting set of partial differential equations defining a pupil mapping system thus 
reduces to: 

(14) 
(15) 
(16) 

(17) 
(18) 



dgi dg2 
dy dx 


= Aii.yy 




dh 


fi{x,y)- 


- X 


dx 


2Z 




dh 


f2{x,y)- 


-y 


dy 


2Z 




dh 


_ x-gi{x 


y) 


dx 


2Z 




dh 


_ y-g2ix, 


y) 


dy 


2Z 





where /i and /2 are the invers e functions associated with gi and g2- Variations of these 
equations were first derived by IVanderbei fc Traubl (120051 ) for circular geometries using ray 



optics. They are the governing relationships for the design of a PIAA unit. 

While it is possible to design pupil remappers with two mirrors of different areas, we 
will restrict our study here to area preserving systems, meaning the set S of (x, y) over the 
entrance pupil is the same as the set of (x, y) over the exit pupils. 
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Note that, as shown in 



Vanderbei &: Traubl (120051 ) . it is sufficient to solve Eqs. [T 



and [12] for gi{x,y), g2{x,y), and h{x,y), since, by conservation of optical path length, 
h{x,y) = h{fi{x,y)j2{x,y)) - ^ + ^ 



(19) 



where /i and /2 are the inverse of gi and g2- This conservation of optical path reduces the 
problem to solving only one second order partial differential equation (by solving Eq. [12] for 
91, Eq. [13] for g2 and substituting into Eq. [Hj). This finishes our derivation of the design 
equations of PIAA systems using ray optics. We now focus on applying them to the case of 
separable apodization functions on square apertures. 



2.2. Mirror shapes for a square geometry 



As we mentioned above, there is no general solution to the Monge- Amp ere equ a tion 



and thus no closed fo r m pu pil mapping design for arbitra ry pupil geometries (iMongel (Il78ll ) 



Tilmann fc Vladimir! fl2002[ )). 



Vanderbei fc Traubl (l2005l ) restricted themselves to a circular 



geometry, which matches a typical telescope entrance pupil. The resulting azimuthal 
symmetry reduces the pupil mapping equations to one-dimensional partial differential 
equation's in radius, leading to relatively simple solutions for the mirror profiles. We 
consider here the case of a square aperture inscribed in the circular telescope pupil. 

To make the design tractable on square pupils we further reduce the problem by forcing 
a separability condition on the prescribed apodisation. 



•-;^2 



A{x,yy = AxixYAyiy) 



(20) 



Under this constraint, there are multiple solutions of Eqs. [12] and [13] We thus choose to 
search for gi, g2 and h such that: 



9i{x,y) 



(21) 
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h{x,y) = hx{x) + hY{y) 



(22) 
(23) 



Note that Eq. [23] can be deduced directly from Eqs. [2T] and [221 These constraints reduce 
Eq. [13] to the separable form, 

dgi dg2 



Ax{x)'Ay{yy 



dx dy 

Considering a pupil of width a, we use the facts that both area and apertures size are 
conserved, namely j Ax{xYdx = a , or j dxi{x) = a, and we find, via separation of 
variables, 



(24) 



9i[x) 

92{y) 

h{x,y) 



Axiafda-"^ 

-a/2 ^ 

y 

a/2 



Avia^da 



a 
2 



Axiafda - u ] diy+ — 



a/2 \J-a/2 



1 rv 
Z 



(25) 
(26) 

Ayiafda - T ) rfr(27) 



a/2 \J -a/2 



We can apply this result to design a PIAA that reproduces an optimal high contrast 
apodisation over square apertures. The idea of using square apertures and taking 
advantage of the highe r cont rast on the diagonal of the PSF was first introduced by 



Aime et al. 



Nisenson &: Papalioliosl (2001 ). who used sonine and cosine apodisations. This idea was 



( I2OOII ) who introduced an interferometric induced apodisation. 



re-visited by 

They also pointed towards u sing optimal functions for square apertures, called prolate 



functions (lAime et al 



( l2002[ )). Here we use a combination of a horizontal and vertical cross 
prolate apodizer. Each apodizing screen yields a 10~^ contrast, resulting in a constrast on 
the diagonal of the PSF of 10~^°. No te that this is also the underlying idea behind the 



design of checkerboard shaped pupils (jVanderbei et al. 



(I2OO4I )). Using prolate apodizations 



in each orthogonal direction, we thus find two mirror shapes. Fig. [3l which lead to a 
PSF that achieves high contrast in a region around each diagonal. Fig. [1] Note that the 
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deformations along the axes of the square apertures are smaller magnitude than in the 
azimuthally symmetric case, which seems natural since separable square apertures do not 
present a full 360 degrees search space. 

2.3. Performance 

Except for the throughput loss associated with inscribing the pupil in a square, a square 
pupil mapper has all of the advantages of classical PIAA coronagraphs. When designed 
using only ray optics it achieves the needed contrast with nearly 100 percent throughput. 
Moreover, because of the unique combination of the close inner working angle of square 
designs and the magnification effects of a pupil mapper system, described in Appendix A, 
it has a very small inner working angle, as close as 1.5A/D. This magnification is due to 
the dependence of optical path length on the angle of incidence. Our pupil mapping design 
achieves a contrast of 10^^° with an inner working angle (IWA) of 1.8A/D. However, like 
other square pupil apodisations, this is not achieved everywhere but only in the region 
around the diagonals covering roughly 63 percent of the image plane. This is illustrated on 
Fig. E] where on-axis and off-axis two dimensional PSFs are shown. 

Unfortunately, like all pupil mapping approaches, this ray optics performance is not 
achievable in practice. Contrast is limited by light diffracted off the edges of Ml. Under 
narrow band illumination the 10"^" contrast can be retrieved using a Deformable Mirror 
based wavefront controller that treats the electrical field oscillations due to the propagated 
edges of Ml as any other chromatic wavefront error. However, because the magnitude 
of such oscillations varies rapidly and periodically with wavelength, it is fundamentally 
difficult for wavefront control technologies to correct them down to high contrast levels, 
eg 10^^°, over a broad spectral bandwidth, AA/A > 0.05. Here we choose to explore a 
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solution, based on apodizing screens, which yields high contrast over a broadband but at a 
cost in throughput. Such designs are intrinsically achromatic and trar isfers the problem o f 



bandwidth to the control of wavelength dependent optical aberration. 



Pueyo fc Kasdin 



Pueyo et al. 



fl2009ah : 



(120071 ). In the next section we present an analytical treatment of the 



edge-diffraction properties of PIAA designs. 



3. Diffraction Analysis 



The previous section derived the pupil remapping equations usi ng only ray optics. A s 



we no ted earlier, it has already been shown for circular geometries (jVanderbei fc Traub 



([20051)) that when a careful diffraction analysis is carried out, the contrast yielded by a pupil 
mapping system is degraded by many orders of magnitude. In this section we present a new 
approximation to the diffraction integral on arbitrary geometries and show how it leads 
to the pupil mapping apodization in the ray optics limit. This approximation is geometry 
independent and has multiple applications beyond the scope of this communication. In 
particular it has been applied to design robust and ef ficient nume r ical p ropagators of 



arbitrary electrical fields through PIAA coronagraphs. 



Krist et al. 



f l2010l ). a tool that is 



critical for the design and diagnostic of experimental validations. 



3.1. The Main Theorem 

The exact electric field at the exit pupil of a PIAA unit can be found using the Huygens 
integral over all the possible optical paths from {x,y) at Ml, to a point {x,y) in the output 
plane, after M2, 

Eoutix,y) = ^ f f e'^Q^-'^'^'^^dxdy (28) 
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where Eq is the constant entrance electric field, Q{x,y,x,y) is the optical path length 
between any two points in the entrance and exit pupils given by Eq. [TJ We have used 
the paraxial approximation to replace the optical path length in the denominator with 
2Z. Unfortunately, this integral is intractable for analysis and computation. The classical 
Fresnel approximation expands the optical path length in the exponent ial about its value 
along the axis of the optical system (x = and y = 0). As was shown iu iVanderbei fc Traub 



( 120051 ) ■ this approximation is too crude for a pupil mapping analysis. Instead, we expand 
the optical path length about the constant Qq along the ray between the pupil remapped 
coordinates, {x,y) and Xi = gi{x,y), yi = g2{x,y) in Eq. [TTl 

Q{x,y,x,y) = Q, + --^^^^^^^^Jx - x.f 



+ 7;l^TL=....=Jy-y^y■ (29) 



2 dy"^ 'x=Xi,y=y, 

Note that the first order derivatives vanished because the functions gi and g2 where chosen 
to minimize the optical path length. This means we can write the Huygens integral in 
Eq. M as 

Eoutii,y) = / / e'^'^'^^^'y^'^yUxdy (30) 



2i\Z 

where G2 only contains terms of the second order, (x — Xj)^, (x — Xi){y — y^), {y — yi)'^. We 
as sume here t hat th e higher order terms are negligible, using the same arguments presented 



m 



Goodman! (119681 ). The field after M2 is thus mainly described by an integral over the 
quadratic terms of the Taylor expansion of the optical path length across the apparatus. 
This brings us to our main theorem: 



Theorem 1 Consider an arbitrary pupil remapping system designed such that the geometric 
PIAA equations are satisfied. Then, as long as the curvature on the mirrors is smooth 
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enough, the propagated field after M2 can be approximated as: 

Eout{x,y) = ——^ 
liXZ 

where Xi = gi{x,y) and yi = g2{x,y) are the inverse mapping at Ml for the independent 
variables {x,y) and Sq = S{gi{x,y), g2ix,y),x,y). 

The proof of this theorem is given in Appendix B. In the remaining sections we will 
use Eq. [31] to compute the diffraction through a square PIAA system. 




ITT ( dfl I 

g ASq V dx I x^ ,yj 



'dxdy (31) 



3.2. Validity of the second order approximation 

A direct consequence of Theorem [T] is a new derivation of Eq. [H]in the ray optics limit. 
If we consider the behavior of Eq. [31] as A — or equivalently D oo, then the asymptotic 
value of the complex Gaussian in Eq. [31] is the determinant of the Jacobian of the change 
of variables defined by the pupil rem apping. Such a result can be rigorously proven using 



the stationary phases approximation (iBorn fc Wolfd (IIQQqI ). We thus find that 



dfidh _ d}2afi W dx dy dy dx 

Y dx dy dx dy 

The goal of the pupil mapping design is to set the field at the exit pupil equal to the desired 
apodization, A{x^y). Equating Eq. [32]to A{x,y) indeed leads to Eq. [S] the Monge-Ampere 
equation. This analytical result not only provides a sanity check that our diffractive optics 
approach asymptotically behaves as ray optics but also provides a unified theoretical 
understanding of pupil mapping systems. 

We now expand our discussion about the validity of this second order approximation at 
visible wavelength and with finite pupil sizes. The functional form of the integral in Eq. [31] 
is similar to the Fresnel approximation: it is composed of a quadratic exponential centered 
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at the location of the incident ray optics wavelets, {xi,yi). The corresponding effective 
propagation distance has been adjusted according to the effective focal length a ssociated 

, 1^. y,. We follow the presentation on pp 66 - 68 of 



with the corresponding ray — 



Goodman 



(119681 ) and adapt it to Eq. |3T1 We are interested in proving that the terms of order higher 
than 2 in the Taylor expansion of the Optical Path Length do not contribute to the value 
of the propagation integral. For a given point (x, y) at M2, these terms scale as: 



27r 

maxMi{^[Q{x,y,x,y) - Qo - G2ix,y,x,y)]} 



TT 



CX 



4Ag; 



x — Xi)+2- 



[x — Xjjiy — yj) + — — 



iy-y^ry} 



dx 

where A^x^^yi) stands for the region of Ml, centered around {xi,yi), that mostly contributes 
to the value of the fiel d at M2. Follow ing the local stationary phase applied to complex 



gaussians described in 



GoodmanI (119681 ). we find that the area of A(xi,yi) is: 
16Ago 





dh 


_ dh 


dfi\ 




^ dx 


dy 


dx 


dy ) 





(33) 



Thus the condition on the higher order terms becomes: 



-27r, 



maxAfi{^[Q(a;,?/,x,?/) - Qo - G'2(x, x, y)]} ~ 167rA < 1 



(34) 



This condition is always true for the PIAA designs considered here. This discussion 
illustrates that the propagation integral through pupil re-mapping aspherical optics 
presented here is, from an analytical point of view, just as accurate as the Fresnel integral 
when propagating through parabolic optics. Since for the remainder of this paper we 
will only focus on analytical evaluations of Eq. [311 discussing its numerical precision and 
performance when integrated using a set of discrete two dimensional arrays is beyond ou r 



scope. This problem has been addressed in a separate commu nication IPueyo et al. 



Moreover we recently published a comprehensive comparison. 



Krist et al. 



the numerical accuracy of the present method and a S-Huygens model. 



(120091 ). 



(2010), between 



Vanderbeil (120061 ): 
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Belikov et al.l ( l2006l ). There we found that using a combination of both methods allowed 
large gains in computation speed and memory allocation at a small cost in numerical 
precision. 



3.3. Numerical propagator for square apertures 

We can now use theorem [T] to compute the diffraction limited PSF for PIAA systems 
over square separable apertures. Because of the separability, the cross term in Eq. ED 
vanishes and we find 

Eout{x,y) = /e^«o^a. U„,A ^) ^y\^,,yi'-y idxdy. (35) 

If we drop the piston term in front of the integral and use Eqs. [21] and to substitute for 
the derivatives of the ray optics remapping we find 

Eout{x,y) = —-^ e^^oAx(-) dx e^^o^ycy) dy. (36) 

Z'lAZ J_D/2 J-D/2 

We then approximate Sq = Z and proceed to the following change of variables in the 
integrals: 



u = \lj^J-T-y(---^(^)) (37) 



^ {y-y^{x)) (38) 



resulting in 



E^Ji, y) = M^^AlM e-"'*'. (39) 

This expression can be easily calculated using the Fresnel Integrals since 

rh 

I e^" du = FresnelC{h) — FresnelC{a) + i {FresnelS{b) — FresnelS{a)). (40) 
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Thus computing Eout{x,y) reduces to functional evaluations of the Fresnel integrals, 
which are tabulated in many computational packages, 



FresnelC(a;) = / cos(— u )du (41) 
Jo 2 

FresnelS(a;) = / sm(—u'^)du. (42) 
Jo 2 

Fig. [6] shows the field after M2 that has been computed using this propagator. 
Diffraction effects result in field oscillations that degrade the contrast. The PSF, for square 
mirrors of size D = 3 cm separated by z = 1 m, is shown in Fig. [TJ where the diffraction 
limited contrast is 10"*^. Note, however, that for a separable square geometry the contrast 
floor due to edge propagation effects, with perfect optics, no wavefront error and no 
apodiser, is already at the 10~^ level, as compared to the 10~^ floor in circular geometries. 
Since we chose to use a separable amplitude profile on the square aperture, with a soft 
target apodisation along each dimension, we can take advantage of the deeper contrast 
property along the diagonal. As a consequence the limit on contrast induced by edge 
propagation effects is lower than for circular apertures. This contrast is sufficient enough for 
the detection of self-luminous Jovian planets, but edge effects need to be mitigated in order 
to enable the detection of exo-planets. Note that for moderate contrast there also exist 
azimuthally symmetric solutions that do yield the required contrast to image self-luminous 
planets in th e infra-red and only re quire one of the apodizers presented in the next section. 



In particular 



Pluzhnik et al 



(120051 1 showed that a system designed according to geometrical 
optics only, but including an apodizer to allow manufacturability of the PIAA optics, can 
be used in conjunction with a single deformable mirror in order to yield a 10'' contrast in a 
20 percent bandwidth. 
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4. Hybrid Apodisers-PIAA Designs on Square Apertures 
4.1. Pre- and post-apodisers 



In order to design a pupil rema pping system that wi 



in the image plane of the telescope, 



Pluzhnik et al 



1 create a broadband 10 null 



( 120051 ) introduced a set of pre- and 



post-apodisers to smooth the edges of the pupil and reduce the diffraction effects. We 
revisit the roles of these apodising screens in light of Theorem [1] These screens have two 
effects: 

• a) Suppress the field oscillations due to diffraction from the edges of Ml. The 
discontinuity at the edges of Ml is softened by slightly oversizing the incoming 
beam and applying a smooth apodisation to the oversized portion of the pupil. This 
soft edge needs to be mapped out of the main beam path in order to preserve the 
high contrast apodisation profile. As a consequence, the outer part of Ml, that is 
oversized by a fraction aui of the pupil diameter, is designed to be parabolic and 
diverging. This spreads the edge oscillations outside the edges of M2, resulting in a 
reduction of the magnitude of the phase oscillations on M2. The functional form of 
the pre-apodiser is completely arbitrary as long as it smoothly goes to zero outside the 
edges of Ml (is continuous with continuous first derivative). In practice, we choose a 
function to facilitate the calculation of the diffraction propagation using the special 
functions FresnelC and FresnelS. Most current PIAA designs use a simple cosine 
apodization. 



1 + ^os[j^^{x + ^i)]) for 2 — - 2 



for _| < X < f (43) 



^ l{^ + oos[ji^^ix-§)]) for 



D (1+^ 



In Appendix C we show how this apodization results in semi-anlaytical solutions for 
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the propagation integral. Current research is being carried out in order to explore 
optimal choices for the pre-apodizer. 

b) Limit the strength of the remapping. By limiting the amount of apodization in 
the ray optics design, the diffractive effects at M2 due to the strong remapping are 
reduced. In practice this is obtained by using a PIAA unit whose net apodisation on 
the edges of the pupil is weaker. This means introducing a saturation at the edges of 
the prescribed apodisation so that its net value never goes below a pre-determined 
value that we call /3. The cancellation of the field at the edges of M2 is then achieved 
using a post apodiser that allows retrieval of the ideal prolate profile Ax{S:). Again, 
how we relax the apodization and the post-apo dizer we use to recove r it is arbitrary. 
Here, we use the simple functions suggested by 



Pluzhnik et al 



((20051), 



Apiaa{x) = r] — — 

\ost{x) = (44) 

where rj is a. normalization parameter that ensures that ApiAAi^Y integrates to unity 
so that area is preseved within the non-oversized portion of the beam. The two free 
parameters, aMi and (3 define the choice of apodising screens. Note that the strip of 
width omi-D at the edge of Ml is remapped into a strip of width aM2-D at M2, and 
clipped out of the optical train, as shown on Fig. [SI 

The relationship between aM2 and r] and the two free parameters aMi and (3 can be 
derived as follows. First, the scaling ri of the prescribed PIAA apodisation is chosen such 
that the area within the actual optical train of the coronagraph is conserved: 

fD/2 

D = ApiAAix^dx (45) 

J-D/2 
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Because we seek a continuous curvature for the surface of Ml we also impose the value of 
the prescribed PIAA apodisation for |x| > D/2 to be such that: 

Ax{Dl2) + f3 [3 
AprAA{x) = ^ ^ (47) 

since Ax{x) is almost zero near to the edge of the pupil. Then, using energy conservation, 
the strip of width aMiD/2 is now remapped at M2 into a strip of width 

o^M2D = = aMi{riY^VD/2. (48) 



This yields. 



"Ml = aM2ApiAA{D/2f = "Af2('7Y^^)^- (49) 



Thus, since /5 << 1 and ?7 — 1, the thin strip at the edge of Ml is remapped into a thick 
area at the edges of M2: umi « o:m2- Note that for a given (3, Eq. HQ] corresponds to 
a one to one mapping between omi and aj\/2- Thus we do not loose generality when, for 
practical reasons, we choose to use aM2 as the free design parameter. 

In the next section we present several square PIAA designs. To do so we first define 
the post-apodiser by choosing aM2 and (3 and then use Eq. |49]to compute aMi and thus 
design the pre-apodiser. 

We derived here the relationships between aM2 and rj and c^mi and (3 under the 
assumption that the pupil mapping optics conserve the area within the apodization region 
seen by the coronagraphs. This implies that the actual physical size of Ml is {1 + aMi)D 
and only its central portion of width D is transmitted to the corononagraph. One might 
want to choose a different convention and use D as the physical size of Ml and scale all the 
quantities accordingly. In that case M2 and the transmissive region of the post- apodiser 
are smaller than D. We summarize both conventions in Table [H In the present paper 
we chose to use the convention in the first line of Table [H Finally, we point out that for 
large pupil sizes, whose edge diffraction effects are relatively small, one can design a PIAA 
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unit with a mi > without using a pre-apodiser to taper the edges of Ml. However in 
an effort to keep the scope of this paper as broad as possible, we present designs that use 
pre-apodisers. 

4.2. Numerical results and performance 

Fig [8] illustrates how the light propagates over a one dimensional aperture for such pre- 
and post-apodised designs. This figure shows true electrical field distributions computed 
using our analytical propagator. For the designs in Fig. [10] and [11] we respectively chose 
«M2 = 4 and /3 = 0.15, which yielded omi = 0.08 and a stroke of 62 /im, and aM2 = 4 and 
(3 = 0.05 which yielded umi = 0.01 and a stroke of 585 fim. The large edge oscillations 
are now remapped very far from the optical axis and thus have little impact on the 
contrast. Combined with the fact that M2 is only reflective for |x| < D/2, this effect 
has a highly beneficial impact on contrast: the edge oscillations are not captured by M2 
at a very small cost in throughput. Because of the trigonometric form we chose for the 
pre-apodiser in this study, predicting the contrast and performances of such PIAA designs 
can be achieved without any numerical integration, using only well established functional 
evaluations. In Figs [TD] and [TT] we show the field at M2 for pre- and post-apodised square 
PIAA coronagraphs with respective pupil sizes of 3 cm and 9 cm. In both cases, the phase 
oscillations at M2 are small enough so that the PSFs exhibit a 10^^'^ contrast on the 
diagonal axis (bottom panels of Figs [TO] and [TTl) . For each geometry, that is, each choice of 
D and z, there is an entire family of apodising parameters (aAfi,/3) that will yield a 10~^° 
contrast. This parameter space will be explored in a future communication. Fig. [T2] and 
Table [2] illustrate three designs we obtained with pupil sizes 9, 3 and 1 cm. 

The contrasts of these design were found at a single wavelength of A = 600 nm. To 
check the broadband performance, we computed the PSFs over a 120 nm bandwidth (a 
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20% band centered on A=600nm) with the same combination of aM2 and /3 previously 
used. For this bandwidth, the 10^^*^ contrast is conserved, as shown on Fig. [121 In general, 
the larger the wavelength considered, the larger the stellar halo is since propagation effects 
scale with D'^/{Xz). A broadband optimization of PIAA will be presented in a subsequent 
communication. 

5. Feasibility of PIAA over square apertures: PIAA with deformable mirrors 

In this section we address some issues regarding the practical implementation of the 
PIAA solutions on square apertures. Since most optical manufacturing tools are built in 
order to accommodate circular apertures, polishing a set of square ground PIAA mirrors 
such as the ones used to generate Fig. [131 with deformations close to 600 fim, is a difficult 
endeavor. Herein we focus on the feasibility of a PIAA design using square optics whose 
optical surfaces can be actuated with very high precision: Deformable Mirrors (DMs). Most 
of the current DMs envisioned for high contrast imaging applications are square, and are 
thus perfectly suited to implement the solutions discussed in this paper. The DM PIAA 
unit presented below was designed under the following constraints: 

• The broadband contrast, 20 percent bandwidth, should be below 10"^''. 

• The peak-to-valley deformation should be below 4 /im for each DM. This number was 
chosen based on the state of the art maximal stroke for MEMS DMs. 

• The mirror curvature should be below the maximal curvature obtained by an actuator 
pushed at maximum stroke. In practice this maximal curvature constraint is already 
achieved by the use of a post-apodiser that enables a saturation of the DM induced 
apodisation. 
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The size of the active region of the DMs is fixed to 1 cm which approximately 
corresponds to the size of the current large stroke MEMS DMs. We discuss here designs 
in agreement with the bottom panel of Fig. [T2] and the bottom line of Table. |21 The free 
parameters are the separation z between the DMs, the profile of the pre- and post-apodisers, 
and the DM deformation. 

We note that the deformation of the mirrors scales as D"^ j z and is a weak function 
of the choice of pre- and post-apodizers. Namely, for given mirror diameter D and mirror 
separation z, and for the functional form of the pre- and post-apodizers in Eq. US] and 
Eq. SH there is a wide range of apodization parameters that will yield deformations below 4 
microns. This means that one can choose a geometry according to the maximal deformation 
constraint, and then proceed by tuning the apodizers so the broadband contrast constraint 
is satisfied. Thus, our first step is to satisfy the sag constraint by selecting an adequate 
geometry, 2; = 1 m in our case. We then proceeded by enforcing the contrast constraint 
using numerical propagators that evaluate the electric field ringing at M2, which is due 
to the edge diffraction between Ml and M2. The magnitude of this ringing is the feature 
that limits the broadband contrast of a given PIAA unit. The purpose of the pre- and 
post-apodizers is to mitigate these chromatic high frequency edge oscillations that are 
remapped near the center of M2. These remapped Fresnel rings scale as \z/ D"^. Since we 
chose \z/D'^ to be about 0.01, the apodizers need to be quite strong in order to mitigate 
for large ringing. Consequently, we chose aggressive apodisers, /3 = 1.3 and = 0.5, 
which limits our throughput to 0.19 and increases our Inner Working Angle (IWA) to 
2.4A/-D. The loss in throughput is mostly due to the fact that DM1 has to be considerably 
oversized, and since /3 > 1, the actual angular magnification differs from the one for the 
ideal ray optics designed system: instead of a value close to 2, the angular magnification 
equals 1.4 along the diagonal. The results are shown Fig. [TH The top of Fig. [H] illustrates 
the deformation of two 1 cm DMs separated by 1 m for the apodisers parameters chosen 
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here. The PSF at three wavelengths across the bandwidth in the visible are shown on the 
bottom part of Fig. [H] . 

Finally we discuss here the practical implementation of the solution proposed in this 
subsection. The actual reflective portion of DM1 has to be of 1.18 cm, with a 1 cm inner 
region that is responsible for the coronagraphic remapping. In that region the maximal 
deformation is 4 /im, and this region ought to be actuated since it corresponds to the 
high optical quality part of the beam that will propagate through the coronagraph. The 
outer edges also exhibit a 4 fim deformation, but with a diverging profile and the optical 
quality of this part of DM1 can be crude since this light will not reach the final science 
camera. DM2 exhibits as well a maximal deformation of 4 fim. Thus we believe that 
the solution presented here is a viable alternative to current PIAA solutions, using soon 
to be available DM technology. This result is of critical importance with respect to the 
feasibility of a PIAA coronagraph. Because any space based coronagraph will require two 
sequential DMs for wavefront correction purposes, we are proposing here to accomplish 
both the coronagraphic suppression and the wavefront correction using the same apparatus. 
This not only simplifies the design of the instrument, but also eliminates the difficult 
problem of polishing static PIAA mirrors to sufficient precision. In the design proposed 
here, since the light remapping surfaces and the actively controllable surfaces are the same, 
the requirements on the static aberrations (deformation of these optics without any DM 
command) are considerably relaxed compared to the case where the remapping optics and 
the DM are separate. 



6. Conclusion 



This paper presents recent advances in our theoretical understanding of the impact on 
contrast of free-space propagation between the two aspherical mirrors designed to generate 
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a Phase Induced Amplitude Apodisation. Our results are two-fold. First, we introduced 
a semi-analytical tool to compute the diffraction limited field after PIAA coronagraphs of 
arbitrary geometries. This methodology has already been applied to a number of numerical 
and experimental studies pertaining to the manufacturability of such coronagraphs and to 
their performances in the presence of wavefront errors compensated by Deformable Mirrors. 
We expect this semi-analytical result to be fully integrated in the design of the broadband 
wavefront control loop in experiments validating the ability of pupil mapping coronagraphs 
to yield the dynamic range necessary to image exo-earths from space at very close angular 
separations. 

The second conclusion of this paper stems from applying this novel diffractive formalism 
to design PIAA units on square separable apertures. In the case of exo-earth imaging, it is 
sensible to explore such geometries since current experimental efforts to integrate together 
coronagraphs and wavefront control are based on square Deformable Mirrors. We illustrated 
the performance of several designs over different pupil sizes, and laid down the theoretical 
foundation for a systematic optimization of such coronagraphs. Finally, we emphasized a 
class of solutions that is particularly interesting: a PIAA unit composed of two Deformable 
Mirrors. Indeed, for small pupil sizes, the amplitude remapping necessary to generate an 
apodisation that produces sufficient contrast for exo-earth imaging can be generated solely 
with DMs. While the throughput of the designs presented here is reduced to 20%, these 
solutions considerably simplify manufacturing and wavefront control of starlight suppression 
systems. We emphasize that the performances of DM based PIAA coronagraphs presented 
in this communication are illustrative realizations of a vast parameter space. Ongoing 
theoretical work is being carried out in order to mitigate the loss in performance associated 
with these solutions. In particular we are currently investigating optimizations that seek 
to maximize throughput, and therefore angular resolution, for given levels of contrast and 
bandwidths. 
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A. Geometrical angular magnification 



We follow the presentation of iMartinache et al.l ( 120061 ) and define the local angular 



magnification at M2 as the spatial frequency of a tilted wavefront at Ml, as seen from M2. 
Assuming that the incoming wavefront is written as: 

Ei„(x, y) = e^lT^Ml(xcos</,+ysin</,) 

Then, because of the pupil remapping, the first order of the outgoing wavefront will be: 

Eout{x,y) = e^ljTAf2(S,!/)(£cos0+ysin0)_ ^^2) 

We define the local angular magnification as 

M{x,y) = ^-^^^ (A3) 
7m 1 

Moreover Eout{x,y) is directly given by the geometrical remapping, 
resulting in, 

M{x,y) = cos 0)2 + (^. sin 0)2 ( A5) 

Here we define the overall angular magnification using the barycenter of the PSF of an 
off-axis source. More realistically it should be defined using the maximum of such a 
PSF, but in practice the two results are very close and the convention chosen in this 
paper allows us to quickly compute A^. Then overall magnification, A4, is defined as the 
average of the local magnification weighted by the pupil transmission function. Carrying 
out the calculation for the design presented associated to Fig [Til for which /3 < 1, yields 
magnifications that oscillate between 2.1 and 2.2 depending on 0, the orientation of the off 
axis source. For simplicity we choose to use the maximum value, Ai = 2.2 
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B. Proof of the analytical propagation theorem 

Our objective is to show that the second order terms in the Taylor expansion of the 
optical path length about the rays of pupil mapping, namely G2{x, y, x, y) in Eq. [301 reduces 
to the exponent shown in Eq. EU For clarity, we begin by repeating this Taylor expansion, 

G2{x,y,x,y) = - x,)^ 

where Xi = gi{x,y) and yi = g2{x,y). Note that in the proof which follows, the coordinates 
{x,y) appearing in Q are independent of {x,y)] the partials are taken only with respect to 
X and y and then evaluated at the pupil mapped locations Xj and yi. 

We proceed by evaluating the three second partial derivatives of the optical path length 
using Eqs. [1] and [2) 



S{x,y,x,y)^ dx^ 



dy'^ 



Qy2 Qy2 



{h 



h) + ^')S{x,y,x,yf-iiy 



S{,x,y,x,yf 



+ 



(B3) 



dxdy dxdy dxdy 

^ if.fy + ih-h)S^)S{x,y,x,yf 
S{x,y,x,y)^ 
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{{x-x) + {h-h)^){{y-y) + {h-h)^) ^ d'h 
S{x,y,x,yy dxdy 

Because of symmetries, we need only carry out the derivation for the second derivative 
with respect to x and the mixed derivative. We re-arranere and as: 

^ ° OX'' axay 

d'Q_l d'h h-h {SfJ'-{{x-x) + f^{h-h)f 



d'Q _ d^h h-h s^fJ^^-{{x-x) + f^{h-hMy-y) + fyih-h)) 

dxdy dxdy^ S ' ^ ' 

Note that we have not yet replaced the coordinates by the actual ray optics values; thus 

Eqs. IB5I and lB6] are valid for all (x, y) and {x,y). 

Next we rewrite Eqs. 15 and 16, which give the height h of Ml as a function of (x, y) 
and the pupil mapped locations fi{x,y) and f2{x,y): 

^ A ^ - h{fi{x,y)j2ix,y)) \ ^ fi{x,y) - x ^^^^ 

dx y S{x,y,fi{x,y),f2{x,y)) J S{x,y, fi{x,y), f2{x,y)) 

^ fl^ K^^y) - Kfiix^y),f2ix,y)) \ ^ f2{x,y)-y ^^^^ 

dy I S{x,y,fi{x,y),f2ix,y)) / S{x,y, fi{x,y)j2ix,y)) 



or, 



dh ~ Bh 

S{x, y, A, f2)^ = -{h{x, y) - M/i, f2))^ -{x- h) (B9) 

dh ~ dh 

S{x, y, A, A)^ = -(/^(x, y) - h{h, A))— -{y- A). (BIO) 



We now use Eqs. IB9I and IBIOI to replace the partial derivatives of h in Eqs. IB5I and 
Beginning with Eq. IB5t we see that it is composed of three additive terms. Evaluating at 
{xi,yi), the third one can be simplified as follows: 

isfj-iix-x) + f^ih-h)y 



(Bll) 



{S{xi,yi, x,y) Q^\xi,yi) (S(xi,yi, x,y) Q^\xi,yi) (u-\o\ 
S{xi,y„x,yy 

(B13) 
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where we have used Eq. IB9I evaluated at Xi = gi{x,y) and i/i = g2{x,y) to cancel the two 
terms in the numerator. In order to simplify the second term of Eq. IB5I we re-write Eq, 
as: 

dh ~ dh 

g^{S{x, y, A, A) + {h{x, y) - h{f„ m = —Q{x, y, A, A) = (A - x). (B14) 
Taking the partial derivative of Eq. IB14I with respect to x yields: 



— 2Q{x^ /i> j2> + - — — = - — 1 



ox'^ ox ox ox 

We now evaluate this equation at Xi = gi{x,y) and yi = g2{x,y). By definition, the ray 



B15) 



optics coordinates are such that ^ 



0. Consequently, 



dx"^ 



, Q{xi,yiJij2) = ^ 

(xi,yi) OX 



or, dividing through by S, 

d'^h, h-h. 



dx"^ 



1 + 



S 



= 1(^ 

{xi,yi) Sq dx 



-1) 



where we define Sn as 



(B16) 



(BIT) 



•^o = S{gi{x,y),g2ix,y),x,y) = S{x„y„ fi{xi,y.i}, f2ixi,yi)). 



(B18) 



Substituting into Eq. IB5I leaves 



dx^ 



So So dx 



-V. 



1 5/i 



So dx 



(B19) 



Turning to Eq. IB61 we see that it is composed of two additive terms. The second one 
can be simplified as follows: 

S'fJ^ -{x-x + f^{h-h)){y-y+ f (/. - h)) 



yi, X, 



dh 
dy 



S^ 



S{xi, yi, X, y)^^ 



Xi,y; 

dh 
dy 



S{xi,yi,x,yY 



(B20) 
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where we have used Eq. IB9I and Eq. IBIOI evaluated at Xi = gi{x,y) and yi = g2{x,y) to 
cancel the two terms in the numerator. In order to simplify the first term of Eq. IB6I we 
again use Eq. IB14I Taking the partial derivative with respect to y yields: 



d'h ..^dhdQ dh 

-K-^Q{x,yJij2) + -K-^ = 
oxoy ox oy oy 



(B21) 



We now evaluate this equation at Xi = gi{x,y) and yi = g2{x,y). By definition, the ray 



optics coordinates are such that ^ 



0. Consequently, after again dividing by 5*, 



d'^h ^ h{x,y) - h{hj2) ^ 
dxdy S{x,y J1J2) 



= 1^ 

^i,yi So dy 



(B22) 



Note that we also could have also started with Eq. IB9I and taken its partial derivative 
as a function of x and found 



dxdy 



Since the condition ^ 

ax 



iS{x,yJij2) + h{x,y) - h{f 1J2)) 



dx '"'"^^ 



(B23) 



dh 
dy 



has to be enforced in order for continuous curvature 



Xi,yi 



solutions for Ml and M2 to exist, these two results are equivalent. 
We now substitute these results into Eq. IB6I to find: 



d^Q 



dxdy 

This finishes our proof since the 

a'^OPL 



1 dh 



So dy 



1 dh 



So dx 



Xi ,yi 



(B24) 



a^oPL 

dy^ 



can be calculated using exactly the same procedure 



as 



C. Semi-analytical propagation of the pre-apodiser through the remapping 

mirrors 

In this appendix we detail how we propagate the pre-apodiser using the special 
functions FresnelC and FresnelS. For clarity we work in one dimension, the full two 



- 33 - 



dimensional result can be obtain using a tensor product in the case of separable designs. 
Our goal is to derive an expression for 

EoutPre{i) = TTT^ / Ep,e(a;)e^^^^^^'"'^^'^ dx (CI) 

that only involves function evaluations of FresnelC and FresnelS. We achieve this by 
choosing a pre-apodiser with the functional form in Eq. 



< X < 



2 ^ - ^ 2 



1 for -f<a;<f (C2) 

^ i(l + cos[^(x-f)]) for f <x<^ 

Using superposition we can write that the field at M2 resulting from the propagation of 
such a pre-apodiser through a pupil mapping unit is the sum of the contribution of the 
right, left, and center portions of this such an apodisation function, 

EoutPre{x) = Eieftix) + Ecenter{x) + ERight{x). (C3) 

We adopt the following shorthanded notation for the Fresnel integrals 

Mx+l/2) 

T,ix) = / e't" du. (C4) 

The contribution of the center portion of the apodiser can be computed easily using the 
results described in this paper 

Ecenterix) = 7 rr u (^). (C5) 

To computed the contribution of the left side of the pre-apodiser we use the Euler formula 

cos[- (x )] = -(e o^An^"^ 2)+g*D.„ii^ 2)y (C6) 

Da Ml 2 2 

For each exponential component we add the exponents of the pre-apodiser and the 
propagation kernel and we complete the squares. We then factor out of the integral the 



- 34 - 



exponential terms that only depend on x and re- write it as a function of J^u{x)- We do not 
detail here the actual derivation and only provide the reader with the final result of such a 
calculation, 

EiefAx) = 7 .Mi^ [—^ + 77 + — ^) (C7) 



i( 2 (-)_±ro=ji^iA^A^ jr.(^^) 1 ^Ml , 2^ Af 1 ) , 

+ °Mi "A/iD ^ 5^ 'jr ^ „ ( ^ ^ -I ^ ^ 



^a^^ ■ 2 ■ 4 ■ 



Similarly we find that 



J^Right[Xj - J- f^ ^MiO [—f^ o ^) y'^^) 

V Az 2Ax(i) JJ Z 4: 

+ e -Ml <^Mii5 '^,/XiW^(^7^ -0-—A n2 ) 

V Az 2Ax(S:) L) Z 4: 1/^ 

As a consequence the propagation of a pre-apodiser that has the functional form prescribed 
by Eq. |43] through the pupil mapping unit can be computed using solely functional 
evaluation of the Fresnel special functions. 
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Off-Axis PIAA Off-Axis PIAA, side view 




On-Axis PIAA, side view On-Axis PIAA, side view 

d=0 d=0 




Fig. 1. — Geometry and notation. Top Left: Three dimensional representation of a pupil- 
to-pupil off-axis PIAA system. Top Right: Side view of the geometrical remapping in a 
pupil-to-pupil off-axis PIAA system. Bottom Left: Side view of the geometrical remapping 
in a pupil-to-pupil on-axis PIAA unit: This is the configuration that is studied in this 
communication. Bottom right: Side view of all the rays contributing to the diffractive field 
at a point of coordinates {x, y) at M2. The ray corresponding to the geometrical remapping, 
which has coordinates {xi{x,y),yi{x,y)) in the input plane, is highlighted. 
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Fig. 2. — Illustration of the pupil mapping effect on the light intensity distribution. Left: 
the energy distribution is uniform on Ml. Right: the energy distribution on M2. The area 
of each square is then proportional to the intensity of the apodisation function A{x, y)"^. 



- 39 - 




Fig. 3. — Apodisation in square apertures. Left: One dimensional apodisation. Right: One 
dimensional shape of Ml and M2. Aperture size D = ?> cm, mirrors distance z = 1 m. 
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Fig. 4. — One dimensional profiles along the horizontal axis and the diagonal of the square 
aperture PSF. The x axis is in units of X/DonSky This PSF was computed assuming that 
ray optics hold perfectly. 
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Fig. 5. — Case of square apertures. Top: Two dimensional PSF. This PSF was computed 
assuming that ray optics hold perfectly. Bottom: A companion of brightness 10"^ was added 
to the top right corner at 2X/D 




Fig. 6. — Propagated field after M2 for a pure PIAA unit. Top: Amplitude. Bottom: Phase. 
The X axis is in units of X/DonSky Aperture size D = 3 cm, mirrors distance z = 1 m 
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Fig. 7. — Diffraction limited PSF. Top: Two dimensional PSF. Bottom: One dimensional 
profiles along the horizontal axis and the diagonal of the PSF. The x axis is in units of 
\/DonSky Aperture size D = ?> cm, mirrors distance z = 1 m 
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Fig. 8. — Electrical field propagated through the PIAA unit. Top panel: transmissivity 
profile of the pre-apodiser. Note that the pupil needs to be oversized by a factor umi in 
order to accommodate for the smooth roll off at it s edges. This results in a loss of throughput. 
Middle panel: effect of the propagation through the PIAA aspherical optics; most of the light 
is concentrated in the center of the beam, the smooth edges of the pre-apodiser help damp 
the edge diffraction ripples. Note that these ripples are spatially extended over large zones of 
dimensions ajv/2-D surrounding the transmissive area of the post-apodiser. Bottom panel: the 
post-apodiser retrieves the expected prolate profile and clips the high frequency edge ripples 
out of the main beam path of the coronagraph. Note that the phase oscillations actually 
transmitted are much smaller that the ones in the outer region that has been suppressed by 
the post-apodiser. While in general 0^2 ^ «mi their values have been exaggerated in this 
figure for illustration purposes 
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Fig. 9. — Mirror profiles in the presence of pre and post apodisers. Tlie top and bottom 
panels are similar to the previous figure. Middle panel: the central beam profile is obtained 
using the area preserving mirror shapes calculated according to the ray optics model and 
using the saturated mirror induced apodisation in Eq. HH The outer portion of Ml is chosen 
to have a constant negative curvature. The focal length of this outer parabolic portion is 
chosen such that the curvature of Ml remains constant. This results in a greatly expanded 
area over which the edge oscillations are located in the plane of M2. The bulk of these edge 
oscillations is then clipped by the post-apodiser. Note that the parabolic outer region of M2 
does not need to be actually manufactured since it does not contribute to the field seen by 
the coronagraph. 




Fig. 10. — Top: Propagated field after M2 for a hybrid PIAA unit. Left: Amplitude. Right: 
Phase. Aperture size D = 3 cm, mirrors distance z = 1 m. Bottom: PSF of a hybrid PIAA 
unit. The throughput of the design shown here is 0.55. 
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Fig. 11. — Top: Propagated field after M2 for a hybrid PI AA unit. Left: Amplitude. Right: 
Phase. Aperture size D = 9 cm, mirrors distance z = 1 m. Bottom: PSF of a hybrid PIAA 
unit. The throughput of the design shown here is 0.89. 
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Fig. 12. — Monochromatic PSF of three different square PIAA designs. The pre-and post 
apodiser have been chosen so that the PSF features a 10^^° contrast along the diagonaL 
Top: 9 cm mirrors, separated by 1 m. Middle: 3 cm mirrors, separated by 1 m. Bottom: 1 
cm mirrors, separated by 1 m. Note that as the pupil size decreases, the apodisers need to 
be stronger. This implies that the throughput decreases and the IWA increases. 
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Fig. 13. — Broadband contrast (540, 600 and 660 nm) for a hybrid PIAA solution with 9 
cm mirrors. This figure illustrates how the design presented on the top panel of Fig. [T2] is 
truely a broadband design 
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Fig. 14. — Top: Shapes of DM1 and DM2 in the case of a PIAA unit with a 1cm pupil. 
Bottom: Broadband contrast for a hybrid PIAA solution with 1 cm mirrors. AA = 120 nm. 
The throughput is 0.19 
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Physical size of Ml 


Size of the useful region of Ml 


Physical size of M2 


Throughput 


{l + aMi)D 


D 


D 


1/(1 + aMi)^ 


D 


(1 - aMi)D 


(1 - OiM2)D 


(1 - OiMlf 



Table 1: Two equivalent conventions to describe hybrid apodised PIAA designs. Note that 
to first order they are equivalent in terms of throughput. In this paper we chose the one in 
the top row 
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Pupil Size 


Maximal Deformation 


Throughput 


Angular magnification 


IWA 


9 cm 


585 /im 


0.89 


2.2 


1.8 


3 cm 


62 nm 


0.55 


2.1 


1.9 


1 cm 


4 iim 


0.19 


1.4 


2.3 



Table 2: Summary of the performances of the designs discussed above. The separation of 
the mirrors is maintained constant at 1 m and the contrast is constrained to be below 10~^° 
along the diagonal within a 20% bandwidth around 600 nm 



